Robust and fragile PT-symmetric phases in a tight-binding chain 
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Abstract 

We study the phase-diagram of a parity and time-reversal (VT) symmetric tight-binding chain 
with TV" sites and hopping energy J, in the presence of two impurities with imaginary potentials 
located at arbitrary (P-symmetric) positions (m, m = N + 1 — m) on the chain where m < N/2. 
We find that except in the two special cases where impurities are either the farthest or the closest, 
the PT-symmetric region - defined as the region in which all energy eigenvalues are real - is 
algebraically fragile. We analytically and numerically obtain the critical impurity potential jpt 
and show that jpt oc 1/N —¥ as N — > oo except in the two special cases. When the P7~ symmetry 
is spontaneously broken, we find that the maximum number of complex eigenvalues is given by 
2m. When the two impurities are the closest, we show that the critical impurity strength jpt in 
the limit N — > oo approaches J (J/2) provided that N is even (odd). For an even N the VT 
symmetry is maximally broken whereas for an odd N, it is sequentially broken. Our results show 
that the phase-diagram of a PT-symmetric tight-binding chain is extremely rich and that, in the 
continuum limit, this model may give rise to new PT-symmetric Hamiltonians. 
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Introduction: Lattice models, including tight-binding chains, have been a cornerstone of 
theoretical explorations due to their analytical and numerical tractability [1], the absence of 
divergences associated with the ultraviolet cutoff (2,3], the availability of exact solutions (4], 
and the ability to capture counter-intuitive physical phenomena including the bound-states 
in repulsive potentials j^] . In recent years, sophisticated optical lattice systems have increas- 
ingly permitted the experimental exploration of lattice models jf], 7|. These lattice models 
are largely based on a Hermitian Hamiltonian. Over the past decade, it has become clear that 
non-Hermitian Hamiltonians with PT-symmetry can have purely real energy spectra j^l, 9] 
and that, when they do, with an appro pria tely re-defined inner product, their eigenvectors 



can be appropriately orthonormalized 



101 ]. The theoretical work has been accompanied, 



most recently, by experiments in optical physics where spontaneous VT symmetry breaking 
in a classical system has been observed in waveguides with VT symmetric complex refractive 
index ll|. In recent years, lattice models with a VT symmetric non-Hermitian "hopping" 



between adjacent levels of a simple-Harmonic oscillator 12j . tridiagonal PT-symmetric mod- 
els 13], and tight-binding models with a Hermitian hopping and PT-symmetric complex 
on-site potential 14j-ll6| have been investigated. These models are of physical interest be- 
cause they lead to unitary scattering even in the absence of a Hermitian Hamiltonian ~v\ . 



In this paper, we analytically and numerically investigate the phase-diagram of a one- 
dimensional tight-binding chain with hopping energy J and two imaginary potentials ±27 as 
a function of the size of the chain N and the positions of the two impurities within the chain 
(m,rh) [if], III]. Our main results are as follows: i) We show that except in the two special 
cases (the impurities are the closest or the farthest), in the thermodynamic limit iV — > 00, 
the critical potential strength vanishes, 7pt(a0 oc 1/N — > where fi = m/N < 0.5 is the 
relative position of the impurity. Thus, for a generic impurity location, the PT-symmetric 
phase in this system is algebraically fragile [18] . ii) The "degree of VT symmetry breaking" , 
defined as the fraction of eigenvalues that become complex, is given by 2fi. iii) When the 
impurities are the closest, the critical potential strength is given by 7pt(/-£ = 0.5) = J 
when N is even and when iV is odd, 7pt(/^ ~ 0.5) — > J/2 as N — > 00. iv) By considering 
the continuum limit of this problem, we argue that the PT-symmetric chain maps onto a 
PT-symmetric continuum Hamiltonian with an imaginary viscous drag term. 



Tight-binding Model: We start with the Hamiltonian for an iV-site tight-binding chain with 
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J > and two impurities with imaginary potentials at sites (m, m = N + 1 — m), 

N 

H = - J ( C n C n+l + C l+l C n) + «7 (<4c m ~ <4 C^J , (1) 

n=l 

where (c n ) is the creation (annihilation) operator at site n, m < iV/2 and m = (iV + 
1 — m) > N/2 is the reflection-counterpart of site m. We only consider the single-particle 
sector and, since periodic-boundary conditions are incompatible with the parity operator, 
our system has open boundary conditions. In the lattice model, the action of the parity 
operator is given by Vc^V = with n = (N + 1 — n), and the action of the anti-linear time- 
reversal operator is given by TiT = —i. The potential term in the Hamiltonian, Eq.fJTJ, 
is odd under individual V and T operations, and hence the Hamiltonian is "PT-symmetric. 
The phase-diagram of this system with the impurities at the end and its equivalent Hermitian 



counterpart have been investigated in detail 



A general single-particle eigenfunction of the Hamiltonian can be written as \tpk) = 
Y^ n =i ^n c li\^) = Y^ n ' l f ) n\ n ) where, according to the Bethe ansatz, the coefficients ipn have 
the form 

Asin(kn) n < m, 

= { Psm(kn) + Q cos(kn) m < n < fh, (2) 
B sin (fen) n > rh. 

Note that open-boundary conditions, along with the requirement H\tpk) = Ek\ipk), constrain 
the eigenfunction coefficients to a sinusoidal form in the regions n < m and n > rh and give 
Efc = —2Jcos(k). By considering the eigenvalue equation at points n = m, m + 1 and their 
reflection-counterparts, we find that the quasimomenta k obey the equation 



M(k) = 



sin 2 (A;(m + 1)) + ( ^ ) sin 2 (A;m) 



sm(k(N - 2m + 1)) 



+ sin 2 (fcm) sin(fc(iV — 2m — 1)) — 2 sin (/cm) sin(/c(m + 1)) sm(k(N — 2m)) = 0. (3) 

The VT symmetry is unbroken provided that Eq.(j3]) has iV real solutions. When 7 = 0, the 
iV distinct solutions are given by the well-known result for a tight-binding chain with open 
boundary conditions, k a = an/(N + 1) where a — 1, . . . , N. Since M(n — k) = (—l) N M(k), 
if ko is a solution of Eq.([3]), then so is (n — ko). It also follows that when is odd, k = tt/2 
is a solution irrespective of the value of 7, and that the corresponding eigenvector has zero 
energy. When m = 1, Eq.()3]) reduces to the criterion for quasimomentum obtained in 



Ref. 15[; in that 7/ J — > 1, the two central k a ~ n/2 become degenerate, the VT 

symmetry is spontaneously broken, and the system develops N — 2 real and two complex 
(conjugate) eigenvalues. 
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FIG. 1. (color online) (a) Left panel shows the permitted values of quasimomenta k/n as a function 
of impurity potential (7/ J) 2 for a chain with N = 20 sites and one impurity at m = 4 (/1 = 0.2, 
thin solid red) or at m = 8 (/i = 0.4, thick dashed blue). The plot is symmetric about k = tt/2. 
As (7/J) 2 is increased, the adjacent quasimomenta k\ ~ ir/(N + 1) and ki ~ 2ir/(N + 1) become 
degenerate (as do their counterparts near k ~ it) and lead to the VT symmetry breaking, (b) 
Right panel shows the critical potential strength 7pp( / u)/J as a function of the chain size N for 
various positions m = fiN of the impurity potential. We see that 7pt(a*) vanishes as N — > 00; 
thus the 'PT-symmetric phase is algebraically fragile except when the impurities are either closest 
to each other or the farthest. 



The left panel in Fig. [T] shows the typical plot of quasimomentum values k^/ir as a 
function of 7, for a chain with N = 20 sites and the first impurity at m=4 (thin solid red) or 
m= 8 (thick dashed blue). Since the plot is symmetric about k = tt/2, we focus only on the 
left-half k/ir < 0.5 and note that since iV is even, there is no solution present at k = 7r/2. 
As the impurity potential 7 is increased, the adjacent quasimomenta and the corresponding 
energy levels become degenerate, leading to the VT symmetry breaking [^.[lol]. We see that 
the critical potential for m = 4 is greater than that for m = 8, 7pt(a* = 0-2) > 7pt(/^ = 0-4) 
and, in contrast with the m = 1 case, the two levels that become degenerate occur in a pair, 
with one near the origin k ~ 7i~/ (N + 1) and other near the zone-boundary k ~ Ntc/(N + 1). 



Therefore when 7(/i) = 7pt(a0 + + , there are four complex eigenvalues. Since the VT 
symmetry breaking can be associated with the development of dissipative channels, we 
define the "degree of PT-symmetry breaking" as the fraction of eigenvalues that become 
complex. For a general m, as 7 is increased beyond 7pt, we find that a total of 2m complex 
eigenvalues emerge and thus the degree of "PT-symmetry breaking is given by 2/z < 1. The 
right panel in Fig. [T]shows the typical evolution of critical potential strength 7pt(a0 with N, 
for jj, = 0.2, 0.25, 0.4, obtained by numerically solving Eq.(j3]). The scaling suggests that the 
critical potential strength for the infinite chain approaches zero, 7pt(a0 oc 1/N — > 0. Thus, 
the VT -symmetric phase, which exists in the region < 7 < 7pp, is algebraically fragile. 
This result can be qualitatively understood as follows: in the limit N ^> 1 and m ^> 1 with 
m/N = /i, Eq.flSJ) can be approximated by 7 2 sin [(1 — 2fi)kN] = 0. If 7 7^ 0, this equation 
has only N(l— 2/x) < iV real solutions, and hence the VT symmetry is broken. We emphasize 
that this argument is invalid when // — > and the corresponding critical potential is given 
by 7pt(a0 — J 15]. It is also invalid when // — > 1/2, the impurities are closest to each other 
and, as we discuss below, the critical 7pt(a0 is nonzero when iV — > 00. Incidentally, we 
point out that the corresponding PT-symmetric phase in a tridiagonal Hamiltonian with 



13|- 



non-Hermitian "hopping" is not algebraically fragile 
Closest Impurities and the Even-Odd Effect: We now consider the case of closest impurities. 
Note that due to the PT-symmetric requirement, when N is even the impurities are nearest 
neighbors with m = N/2, whereas when iV is odd the impurities are next-nearest-neighbors 
with m = (N — l)/2. We will first focus on the case with an even N. In this case, the 



condition M(k) = from Eq. (EJ) reduces to the following equation 

J 2 sin 



(J 2 -7 2 ) S in 2 (^V (4) 



/. I V • I , , 

2 J\ v ' V 2 J 

When 7 = Eq.(jH) has distinct solutions given by k a = a7c/(N + 1). As 7 increases 
the adjacent k a approach each other and when 7 = J, Eq. (J5J) has A^/2 doubly-degenerate 
solutions given by k n = 2nir/(N + 2) where n = 1, . . . , (N/2). When 7 > J, it is clear 
that Eq.(p0) has no real solutions. Thus the VT symmetry is maximally broken and all N 
eigenvalues simultaneously become complex. When is odd, the impurities are at sites 
N ± 1 where Ao = (N + l)/2 is the site at the center of the chain. The equation M(k) = 
then reduces to 



cos(/c) 



sm\kN )+ [^) sin 2 (k(N - I)) 



sm(fcAo)sm(fc(Ao - 1)). (5) 
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This equation has all real solutions provided 7/ J < 1/(2 cos kd) where tt/(N + 1) < kd < 
27r/(N + 1) corresponds to the first degenerate quasimomentum. Therefore, we find that 
as 7 is increased from zero the "PT-symmetry breaks at 7pp = J/2 in the limit iV — > 00 
when adjacent k a near the origin (and their counterparts near the zone boundary) become 
degenerate. Hence, for 7 = 7pp + + , there are four complex eigenvalues. On the other 
hand, Eq.Q has only one real solution, k = tt/2, when 7/J > 1/(2cos&d) where ko < tt/2 
is the degenerate quasimomentum closest to the zone center k = n/2. Hence, the number 
of complex eigenvalues increases monotonically and when 7 > J/ (2 cos ku) ~ J(N + l)/37r, 
it saturates to 2m = N — 1 [if]]. Figure [2] shows the quasimomenta k a (j) for a chain with 
N = 13 and N = 21 lattice sites obtained from Eq.(jS]), and confirms these results. Since the 
VT symmetric nature of the potential dictates the minimum distance between the impurities 
when N is odd or even, the phase-diagram of the chain is sensitive to it even as iV — > 00. 




FIG. 2. (color online) Permitted quasimomenta k(j)/7r for a chain with TV = 21 (thick dotted red) 
and iV = 13 (thin solid blue) sites as a function of impurity strength (7/J) 2 when the impurities 
are closest to each other. As 7 is increased, adjacent quasimomenta become degenerate leading to 
a spontaneous 'PT-symmetry breaking. As N — > 00, we find that the critical potential strength 
7pt -> J/2. 



Numerical Results: We start this section with results for the critical potential strength 
7pr(/ i ) as a function of the relative impurity site location < \x = m/N < 1/2 obtained 
by numerical diagonalization of the Hamiltonian Eq.([T]) for various lattice sizes N, even 
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and odd. The left panel in Fig. [3] shows that, for an even N, apart from finite-size effects 
that are prominent near ft = 1/4 and are also present in solutions to Eq. fl3]), the critical 
potential strength 7pp(/i) is vanishingly small except at fx = 1/N and fi = 0.5. In both 
special cases jpt = J- The right panel in the same figure shows results for odd N. When 



fi = 1/N or equivalently m = 1, we recover the result 7pt = Jyl + l/N [15l|. As in the 
case with even N, we find that 7pt(/-0 is suppressed with increasing N everywhere except 
when fi = 0.5 — 1/2N or equivalently m = (N — l)/2. These results are consistent with 
those obtained through the analytical treatment earlier. 



Tpt/J 




7pt/ j 




FIG. 3. (color online) (a) Left panel shows the critical potential strength ^pt/ J as a function of the 
location fi = m /N < 1/2 of the first impurity obtained via numerical diagonalization for a chain 
with even N. Except at the end-points, fi = 1/N and ft = 0.5, as N increases 7pt(m) decreases. At 
the end points, we find that 7pp = J is independent of the value of N. (b) The right panel presents 
similar results for an odd N, showing that the critical 7pr(^) ck 1/N for all values of fi except for 
fi = 1/N and fi = 0.5 - 1/2N. When fi = 1/N, we recover the result j PT = Jy/l + 1/N Q, 
and when the impurities are next-nearest-neighbors, we find that the critical potential strength 
7pt ->• J/2 as iV -> oo. 



We now briefly explore the change in a (typical) eigenfunction ipk{ n ) as a function of im- 
purity potential 7 in the case of nearest-neighbor impurities (even N). In the PT-symmetric 
phase, an eigenfunction is given by ip<(n) = Asin(kn) for n < N/2 and ip>(n) = B sm(kn) 
for n > N/2 where A; is a quasimomentum that satisfies Eq. (l4"]) . Using the eigenfunction 
constraints and Eq.(j4]), it follows that 

"sin(A;(l + JV/2)) 



B = A 



sin(kN/2) 



, - 7 



A exp (id. 



1) ' 



(6) 
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where the angle # 7 satisfies tan# 7 = 7 sin(A;A^/2)/Jsin [k(l + N/2)]. Figure H] shows the 
amplitude \ipk{n)\ and the phase $(n) of the ground-state wavefunction of a chain with 
N = 20 sites and nearest-neighbor impurities. The top (blue) panel shows that when 
7 = 7ft = J, the wavefunction amplitude is even about the center of the chain and the 
phase is given by # 7 = tt/2, as is expected from Eq.([S]). The bottom (red) panel shows 
that when 7 = 1.01J > Jpt/J, the broken VT- symmetry is reflected in the asymmetrical 
wavefunction amplitudes and in the position-dependent phase factor $(n). These are generic 
features of the broken PT-symmetry phase. We also note that in the continuum limit, the 
eigenfunction ipk{x) becomes discontinuous at the center of the chain while the probability 
amplitude remains continuous. 

Conclusion: We have investigated the phase diagram of an iV-site one- dimensional chain with 
a pair of complex PT-symmetric impurities located at sites (m, m) within it. A remarkab 



feature of such a Hamiltonian, in contrast to a tridiagonal Hamiltonian with real entries [12], 
is that in the "PT-symmetric region, its spectrum remains confined within the energy band 
±2J of the model in the absence of impurities; as the impurity potential 7 is increased, 
the level spacing between adjacent energy levels decreases. Our results show that the VT- 
symmetric phase of such a chain is algebraically fragile except when the impurities are 
farthest from each other or are closest to each other. In the latter case, we find that the 
"PT-symmetric phase survives when 7 < 7pt = ^ (even N) or 7 < 7p^ = J/2 (odd 
N). We note that such a chain offers tremendous tunability due to its variable critical 
impurity strength 7pt(/-0 fo r a finite N, and the corresponding variable fraction of complex 
eigenvalues 2/i which translates into the number of dissipative channels in both classical [ll| 
and quantum systems. Thus, a physical realization of such a model [l6| may offer the ability 
to engineer the level-spacings and the dissipation in this system. 

We conclude by briefly pointing out the continuum limit of this system. In the continuum 
limit, the lattice spacing a vanishes and the number of lattice sites N diverges such that 
the length of the chain L = Na remains constant. Note that since X^C^ptVOO) * s a con ~ 
stant, where Vpt is the PT-symmetric potential, the corresponding continuum Schrodinger 
eigenvalue equation is given by 

d 2 d 
- ~ iF5(x)—'ip k (x) = k 2 ip k (x), (7) 

where the dimensionless impurity strength r = 7/ J for nearest-neighbor impurities (even 
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FIG. 4. (color online) The top (blue) panel shows the amplitude | ipk( n )\ (left) and the phase 
$(n)/7T (right) of the ground-state wavefunction of a PT-symmetric chain with N = 20 sites and 
nearest-neighbor impurities with strength 7 = ^pt = J ■ As expected of a PT-symmetric state, 
the amplitude is even around the center of the chain, and the effect of a non-zero 7 is manifest in 
the discontinuous change in the phase, with $(n) = for n < N/2 and <3?(re) = ir/2 for n > N/2, 
consistent with Eq.([6|). The bottom (red) panel shows the same state when 7/ J = 1.01 and the 
PT-symmetry is spontaneously broken. The broken symmetry is manifest in the asymmetrical 
wavefunction amplitude (left) and a position-dependent phase <E>(n). 



N), T = 27/ J for next-nearest-neighbor impurities (odd N), and the eigenf unctions obey 
boundary conditions ipk{% — ±L/2) = 0. Note that the continuum Hamiltonian is PT- 
symmetric, but not Hermitian, due to the non-commuting parts, 5(x) and —id/dx, of the 
"potential" term. Our results imply that the PT-symmetric phase of this Hamiltonian 
survives as long as T < 1, irrespective of whether the number of sites N in the chain is 
odd or even. Indeed, Eq.([7|) suggests a new class of PT-symmetric Hamiltonians with a 
"viscous drag potential" term of the form V = —ihf(x)d/dx that is not Hermitian but is 
PT-symmetric provided f(x) is an even function of x. In the lattice model, such a potential 
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will correspond to PT-symmetric impurities at multiple locations. Detailed investigation 
of such models will improve our understanding of "PT-symmetry breaking in discrete and 
(classical) continuum systems, that can be realized in optical lattices and waveguides with 



complex refractive index, respectively 
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